Ultimate Turbulent Taylor-Couette Flow 
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The flow structure of strongly turbulent Taylor-Couette flow with Reynolds numbers up to Ra = 
2 ■ 10 6 of the inner cylinder is experimentally examined with high-speed particle image velocimetry 
(PIV). The wind Reynolds numbers Re w of the turbulent Taylor- vortex flow is found to scale as 
Re w cx Ta}/ 2 , exactly as predicted [lj for the ultimate turbulence regime, in which the boundary 
layers are turbulent. The dimensionless angular velocity flux has an effective scaling of Nu u cx 
To ' 38 , also in correspondence with turbulence in the ultimate regime. The scaling of JVti u is 
confirmed by local angular velocity flux measurements extracted from high-speed PIV measurements: 
though the flux shows huge fluctuations, its spatial and temporal average nicely agrees with the result 
from the global torque measurements. 
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The Taylor-Couette (TC) system is one of the funda- 
mental geometries conceived in order to test theories in 
fluid dynamics. Fluid is confined between two coaxial, 
differentially rotating cylinders. The system has been 
used to measure viscosity, study hydrodynamic instabil- 
ities, pattern formation, and the flow was found to have 
a very rich phase diagram Q. In the fully turbulent 
regime, the focus up to now has been on global transport 
quantities 0-0|, which can be connected to the torque 
r, which is necessary to keep the inner cylinder rotating 
at constant angular velocity. In ref. Q the analogy be- 
tween the angular velocity flux in TC turbulence and the 
heat flux in Rayleigh-Benard (RB, see ref. flow was 
worked out, suggesting to express the former in terms of 
the Nusselt number, Nu^, which in ref. Q was found to 
have an effective scaling Nu^ cx Ta 038 with the Taylor 
number (the analog to the Rayleigh number Ra in RB 
flow). Such effective scaling Nu cx Ra 39, characterizes 
the so-called ultimate scaling regime in RB flow |10h12|. 
Following these papers, Grossmann and Lohse |l[ have 
interpreted this scaling as signature of turbulent bound- 
ary layers. They derived Nu cx Ra 1 / 2 x log-corrections 
(RB) and N UuJ cx Ta 1 / 2 x log-corrections (TC). The log- 
corrections imply the effective scaling law exponent of 
0.38. They also made a prediction for the accompanying 
scaling of the wind Reynolds number Re w , namely 



Re w cx Ra}l 2 and Re w cx Ta 1 ^ 2 



(1) 



for RB and TC turbulence, respectively. Here the loga- 
rithmic corrections remarkably cancel out, in contrast to 
what Kraichnan had predicted earlier, namely 



Re w cx Ra 1/2 (log Ra)- 1/2 or Re w cx Ta 1/2 (log Ta)- 1/2 , 

(2) 

which leads to an effective scaling exponent of about 0.47 
in the relevant turbulent regime. In order to verify the 
interpretation of ref. [H and to check the prediction (JlJ, 
local flow measurements are required to extract the wind 
Reynolds number Re w . However, what happens locally, 



inside the TC flow, has up to now only been studied 
for relatively low Reynolds numbers Re < 10 5 , and has 
been restricted to flow profiles and single-point statistics 

In this paper we supply local flow measurements from 
high-speed particle image velocimetry (PIV) at strongly 
turbulent TC flow. From these we will verify that indeed 
Re w cx Ta 1 ! 2 . In addition, from the PIV measurements 
we are able to also extract local angular velocity fluxes. 
These are found to strongly fluctuate in time, but when 
averaged azimuthally, radially, and in time, for the lower 
Ta show a slight axial dependence, which we interpret as 
reminiscence of the turbulent Taylor vortices, and which 
nearly vanishes for the largest Ta we achieve. 

The apparatus used for the experiments has an inner 
cylinder with a radius of r, = 0.200 m, a transparent 
outer cylinder with an inner-radius of r a = 0.279 m, re- 
sulting in a gap- width of d = r a — r.; = 0.079 m and a ra- 
dius ratio 1] = ri/r Q — 0.716. The height is L = 0.927 m 
implying an aspect ratio of T = L/(r — n) = 11.7. More 
details regarding the experimental facility can be found 
in ref. [24(. Here we focus on the case of inner cylin- 
der rotation and fixed outer cylinder. The local veloc- 
ity is measured using PIV. Wc utilize the viewing ports 
in the top plate of the apparatus to look at the flow 
from the top. The flow is illuminated from the side us- 
ing a pulsed Nd-YLF laser (2fl|, creating a horizontal 
laser sheet. The working fluid (water) is seeded with 
20 urn polyamide seeding particles, and is recorded using 
a high speed camera [26j. The PIV system is operated 
in double-frame mode which allows us to have a At far 
smaller than 1//, where / is the frame rate. The PIV 
measurements give us direct access to both the angular 
velocity oj{6, r, z,t) = ug(6,r,z,t)/r and the radial veloc- 
ity u r (6,r, z,t), simultaneously. 

From the latter we extract the wind Reynolds number 
as Re w := uf A d/v, where uf d is the standard deviation 
of the radial velocity. In fig. [l]i?e Ii , is shown as a function 
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FIG. 1: (color online) Re w vs. Ta. The data from re- 
peated experiments at mid-height are plotted as separate 
(blue) dots, showing the quality of the reproducibility and 
the statistical stationarity of the measurements. We have av- 
eraged azimuthally, over time, and in the bulk flow (0.23m < 
r < 0.25m). The straight line is the best fit Re w — 
0.0424ra°- 495±0 010 and the (red) dashed line is the Kraich- 
nan prediction []]§] eq. ((2J). The inset shows the compensated 
plot Rew/Ta 1 ^ 2 vs Ta. The horizontal (green) line is the 
prediction {TJ of ref. Q}. 



of the Taylor number 

Ta = \v{r - ri) 2 {r % + r ) 2 (uji - u ) 2 jv 2 . (3) 

In refs. [1, H| Ta had been suggested as most appropriate 
independent variable of the TC system in order to work 
out the analogy with RB. Here a = ((1 + rj) / '(2^/77)) 4 can 
be interpreted as geometric "Prandtl number" [8(, oj^ is 
the angular velocity of the inner and outer cylinder, re- 
spectively, and v is the kinematic viscosity. Note that 
Ta oc (uji — u! ) 2 : while Ra in RB convection is propor- 
tional to the temperature difference times the given grav- 
ity force, Ta in TC flow is proportional to the angular ve- 
locity difference Wi — lu times the centrifugal force, which 
itself is also proportional to uii — w , implying the square- 
dependence. Therefore, by definition, the two control 
parameters i?e.; (refering to the imposed azimuthal ve- 
locity) and Ta are connected by Rei ~ Ta 1 / 2 , but such a 
trivial relation of course does not exist between the wind 
Reynolds number Re w and Ta (which is a response of 
the systems and refers to the radial velocity) . 

Fig. [T] reveals a clear scaling of the wind Reynolds 
number with the Taylor number, namely Re w oc 
Ta°' 495±0 01 °, which is consistent with the prediction [l| 
Re w oc Ta 1 / 2 for the ultimate TC regime, but incon- 
sistent with Kraichnan's earlier prediction (|2|) of a scal- 
ing exponent 1/2 with logarithmic corrections [l3j]. For 
comparision, we included this relation into fig. [TJ which 
clearly is inconsistent with the experimental data. We 
stress that the cancellation of the log-correction for Re w 
as suggested in [l[ is highly non-trivial and that in RB 



flow in the non-ultimate regimes the wind Reynolds num- 
ber scales as Re w ~ Ra 0A4 (2?]], pronouncedly differ- 
ent than the 1/2 exponent we find here in the ultimate 
regime. Only very recently the wind Reynolds number 
scaling in ultimate RB flow could be measured, also find- 
ing Re w ~ Ra 1 / 2 [HI as predicted in ref. fl. 

Next, as the PIV measurements give us both the 
angular velocity u}(9,r, z,t) and the radial velocity 
u r (9,r, z,t), we can directly calculate the (total) angu- 
lar velocity flux (convective + molecular) 

J"(6,r,z,t) := r 3 (u r u - vd r u) , (4) 

which is made dimcnsionlcss with its value for the lam- 
inar infinite aspect ratio case, J£ m = 2vr 2 r 2 (ui — 
cj D )/(r 2 — rf), giving [H the local "Nusselt number" 

Nu u {6,r,z,t) = J"{6,r,z,t)/J^. (5) 

Indeed, as shown in ref. @, the angular velocity is 
the relevant quantity transported from the inner to 
the outer cylinder, as its flux (|4]) is radially conserved, 
once it is averaged azimuthally, axially, and over time, 
4- ( J w (#, r, z, t)) e z t = 0. In the turbulent regime the 
convective term is the major contributor to the flux in 
the bulk 

In fig. [2^, we show a snapshot of Nu(9, r) at mid-height 
z = L/2 for Ta = 1.5 ■ 10 12 . The quantity shows huge 
fluctuations, ranging from +10 5 to — 10 5 and beyond, 
whereas the average (Nu UJ (9, r 7 t)) e rt = 325 is very close 
to the value Nu^! ob = 326±6 obtained from global torque 
measurements [6| . The local flux can thus be more than 
±300 times as large as the mean flux. Large fluctuations 
have also been reported for the local heat-flux in RB flow 
[30^ . but in that case the largest fluctuations were only 
25 times larger than the mean flux. 

After azimuthal and time averages, (Nu UJ (9,r,t)) g t , 
the fluctuations nearly vanish, see fig. [2]b (revealing some 
radial and height dependence for fixed Ta = 1.5 • 10 12 , 
presumably reminisent of the Taylor vortices) and fig. [5J:, 
where we show the local angular velocity flux r'-profiles 
for rotation rates from u!i/(2ir) = 0.5 Hz to 20 Hz, cor- 
responding to Ta = 3.8 • 10 9 to 6.2 ■ 10 12 . Each profile is 
based on azimuthal averaging, radial binning, and aver- 
aging over 3200 frames (corresponding to 25.6 rotations 
for the three lowest rotation rates, and 32, 64, and 128 
rotations for the fastest rotations rates). For each ro- 
tation rate repeated experiments have been performed 
and the profiles are reproducible. Only in one case the 
turbulent Taylor vortex flow seems to be in a different 
state(s). From fig. we conclude that the spread in the 
repeated experiments decreases with increasing Ta, for 
which the Taylor vortex structure will be more and more 
washed out. In addition, for increasing Ta, not only do 
we measure during more revolutions, but also the trans- 
verse velocity increases, both improving the statistics. 
The dashed lines in fig. [^t correspond to the measured 
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FIG. 2: (color online) (a) Snapshot of the instantaneous convective angular velocity flux, measured at z — L/2, for Ta = 
f .5 ■ 10 12 . The (r, #)-plane and time averaged flux is found to be equal to {Nuu) g , = 325. A corresponding movie is available 
as supplementary material, (b) Local normalized convective tj-flux as functions of r' = (r — r;)/(r — n) for 6 heights varying 
between 0.5 < z/L < 0.73, for Ta = 1.5 ■ 10 12 . The black solid line is the average of the 12 experiments {(Nu^e^^j), which is 
very close to the expected value 326 ±6 from global torque measurements [f|. (c) Local normalized convective angular velocity 
flux vs. radial position r' for various rotation rates, measured at z — L/2. From bottom to top we have repeated experiments 
for Ta = 3.8 • 10 9 , 1.5 ■ 10 10 , 6.2 • 10 10 , 3.8 ■ 10 11 , 1.5 • 10 12 ^,nd 6.2 ■ 10 12 . The dashed lines for the three highest Ta represent 
the Nu u value derived from global torque measurements [g. 



global transport for the three highest rotation rates; these 
values were obtained from the torque measurements @ 
and show already good agreement with our local mea- 
surements. 

An additional axial average is necessary to obtain the 
exact relation between Nu u and the global torque r re- 
quired to drive the inner cylinder at constant velocity 

T = 2TLpJ% m (Nu u ) gtZtt . (6) 

It is the lack of sufficient axial averaging, which accounts 
for the small deviations between (Nuu(6, z,t)) s r t and 
Nu^ ob . Indeed, due to the Taylor-vortex structure of 
the TC flow one would expect some axial dependence of 
{Nuu{0, r, z, t)) g r t , which should become weaker with 
increasing degree of turbulence and thus increasing Ta, 
just as fig-dt suggests. This picture is confirmed in figure 
13 Here we present local measurements of the convective 
angular velocity flux for varying rotation rates, resulting 
in a Taylor number range of 3.8 • 10 9 - 6.2 • 10 12 . For 
each Taylor number we performed multiple experiments 
and measured the Nu u transport at mid-height. The 
blue points are results obtained from PIV measurements 
at mid- height, where the length of the bars indicate the 
error obtained from the repeated experiments. The green 
and orange points are repeated measurements at z = 
L/2 + d/2 and z = L/2 + d, respectively. An effective 
scaling Nu u cx Ta 0A5±0M is revealed for the blue data 
points, while a scaling of Nu u oc Ta 39±0 08 is revealed 
for the orange data points. 

It is remarkable how the flow provides angular velocity 
transport from the inner to the outer cylinder, in spite of 




Ta 

FIG. 3: (color online) Local convective angular velocity flux 
as a function of Taylor number. The blue dots are re- 
sults obtained from PIV measurements and show a scaling of 
Nu„ oc Ta 0A5±0M . The green and orange dots are repeated 
measurements at a height of z = L/2 + d/2 and z — L/2 + d, 
respectively. The black data points are obtained from global 
torque measurements and show a scaling that is less steep: 
Nu u oc Ta 38 . The dashed green line is obtained by match- 
ing two log-layers 0], and has a slope of 0.37 at Ta = 10 9 , and 
0.41 at Ta = 10 13 . The red line is from the turbulent bound- 
ary layer theory of ref. [l]|. It has a slope of 0.43 around 
Ta = 10 9 and 0.44 around Ta = 10 13 . Dark red data points 
are obtained by means of global torque measurements [f| . 



the fluctuative nature, which are seen in figure [2^,. In fig. 
@] we provide a statistical analysis of these fluctuations: 
While the probability distribution functions (PDFs) of 
the angular velocity (fig. QJi) and the radial velocity (fig. 
UJd) are nearly symmetric, the PDF of their product 
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FIG. 4: (color online) Results of three experiment with vary- 
ing rotation rate resulting in Ta = 3.8 -10 11 , 1.5 -10 12 , and 
6.2 • 10 12 , colored in red, green, and blue, respectively. All the 
data shown is averaged over the region 0.23 m < r < 0.25 m, 
and measured at mid-height. All quantities with tildes are 
standardized (shifted and scaled such as to have zero mean 
and unit variance), (a) PDF of the standardized angular ve- 
locity, (b) PDF of the standardized radial velocity, (c) Stan- 
dardized normalized local convective angular velocity flux 
PDF. (d) Cross-correlation coefficient of the angular velocity 
and the radial velocity, the dimensionless decaying time (in 
number of rotations) is found to be 0.07. The corresponding 
length scale can be found by multiplying this number with 
the circumference of the inner cylinder giving 5 = 88 mm, 
which is of the same order of magnitude as the gap width 
d — 80 mm. 



r 3 u r ui oc Nuu (tig. Hfc) is clearly positively skewed. In- 
deed, the cross-correlation coefficient of u r and io (fig. 
|4ji) is relatively large. 

We note that thanks to the PIV measurements of the 
full velocity field, the extraction of the local angular ve- 
locity flux Nu u (9,r, z,t) oc uu r is easier in TC as com- 
pared to the analog temperature flux Nu(x, t) oc Tu z 
in RB flow: in order to obtain this latter quantity lo- 
cally, one has to measure the temperature T and the 
velocity simultaneously. Because a high-precision field- 
measurement of the temperature is presently not possible 
and thus not available, the best one can do for RB flow 
is to measure Nu(x, i) p oint by point (30L 31 1 or use an 
instrumented tracer [32§. 

In conclusion, from high-speed PIV measurements we 
have found the wind Reynolds number in strongly tur- 
bulent TC flow to scale as Re w oc Ta a495±0010 , in ac- 
cordance with the theory of ref. [l| and in conflict with 
Kraichnan's [l3| prediction (|2|) . In addition, we extracted 
the local angular velocity flux and found that Nu^ oc Ta 1 
with 7 ~ 0.39 — 0.45 depending on the axial position 
and consistent with earlier global torque measurements 
@, 0j- For increasing Ta, a small axial dependence of 
Nu u is fading away, reflecting the decreasing importance 



of the Taylor vortices. The next step will be to provide 
full velocity and angular velocity profile measurements, 
including those in the boundary layers, and to extend 
the present measurements to the counter-rotating case 
and other radii ratios ?y, in order to further theoretically 
understand the local flow organization and the interplay 
between bulk and boundary layers in turbulent TC flow. 
A further highly interesting support for the presented 
idea of the close correspondence between the TC angular 
velocity transport in the studied Ta-range with the ulti- 
mate range of RB thermal convection is to identify the 
onset of this ultimate range when increasing Ta; here we 
expect a change of the Nu^ scaling exponent and also a 
transitional change in the widths and profiles of the BLs. 

This study was financially supported by the Technol- 
ogy Foundation STW of The Netherlands. 
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